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Abstract 

The notion of semi-classical states is first sharpened by clarifying two issues that appear to 
have been overlooked in the literature. Systems with linear and quadratic constraints are then 
considered and the group averaging procedure is applied to kinematical coherent states to obtain 
physical semi-classical states. In the specific examples considered, the technique turns out to 
be surprisingly efficient, suggesting that it may well be possible to use kinematical structures to 
analyze the semi-classical behavior of physical states of an interesting class of constrained systems. 
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I. INTRODUCTION 



In the standard Hamiltonian descriptions of fundamental interactions, the canonical vari- 
ables are subject to constraints. Notable examples of systems of this type are gauge theories, 
general relativity and supergravity. In the gravitational case, a key question faced by any 
background-independent approach, such as loop quantum gravity, is whether specific con- 
structions used to impose constraints in the quantization procedure lead to a theory with 
'a sufficient number of semi-classical states'. To analyze this issue one needs a framework 
which spells out the precise meaning of the term 'semi-classical states', introduces strategies 
to construct them and provides tools to analyze their properties. The purpose of this paper 
is to propose such a framework and illustrate its use with simple examples. 

We will restrict ourselves to systems in which the kinematic phase space T — i.e., the 
initial phase space on which constraints are ignored — is a linear space. However, our 
considerations can be extended to cases in which T is a convex subset of a vector space or 
an affine space, and these cases cover a very large class of systems of interest. Our second 
restriction is more significant: We will only consider those constraint surfaces l\i n which 
are level surfaces of linear or quadratic functions on T. This assumption will enable us to 
perform explicit calculations. The Gauss constraint of the Maxwell theory, the constraints 
of linearized gravity, the diffeomorphism constraint of full geometrodynamics, as well as a 
number of constraints often used in the literature for systems with a finite number of degrees 
of freedom are of this form. However, our restriction excludes several important cases, most 
notably the Hamiltonian constraint of general relativity. 

The basic ideas can be summarized as follows. We will begin with a set of Dirac observ- 
ables Oi on T which is sufficiently large to separate points of the reduced phase space T. We 
will then choose coherent states \l/ 7 in the kinematical Hilbert space 7ikm which are peaked 
at points 7 of T in the following sense: Expectation values of the quantum observables Oi 
are the values Oi{^f) assumed by the classical observables Oi at 7, with suitably small fluc- 
tuations. The idea is to restrict 7 to lie on the constraint surface l\m and construct physical 
semi-classical states \1/P hy by averaging \I/ 7 over the the group generated by the quantum 

constraint operator C on 7ikim calculate the expectation values and fluctuations of Oi in 
\1/P hy , and compare them with those in \I/ 7 . 

While the strategy seems natural — even obvious — at first, a priori it is not clear that it 
would be useful. First, the group averaging procedure p need not result in a state which 
has finite and positive norm. Second, even if it does, and therefore defines a state 
in the physical Hilbert space 7i p h y , this state may not at all be semi-classical. For, the 
original coherent state \l/ 7 in the kinematic Hilbert space is not infinitely peaked on the 
constraint surface f. Indeed, the coherent state \l/ 7 will necessarily have fluctuations about 
the constraint surface which, upon group averaging, could contaminate the expectation 
values of the Dirac observables and enhance their fluctuations uncontrollably. Surprisingly, 
this does not happen, at least in simple examples, even when the constraint is quadratic 
and the reduced phase space is a genuine manifold so that H p h y does not admit any direct 
analogs of the kinematic coherent states \l/ 7 . These examples suggest that there may well 
be a general and interesting structure governing the interplay between semi-classicality and 
group averaging and a more extensive and systematic study of the issue may be worthwhile. 

The paper is organized as follows. Section |H] introduces the notation and briefly recalls 
a few properties of the standard coherent states that are repeatedly used in the rest of the 
paper. Section UTTI points out certain subtleties associated with the notion of semi-classicality 
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that are generally ignored in standard treatments but which are often important in practice, 
particularly when dealing with constrained systems. Linear constraints are discussed in 
section IIVI In this case, the expectation values and fluctuations of a very large class of 
Dirac observables in the kinematic coherent states and their group averaged versions 
,j,phy are identical An interesting element is that this result holds even if the fluctuations of 
the constraint operator in the state \l/ 7 are arbitrarily large and indeed even when the point 7 
does not lie on the constraint surface F. In section we consider quadratic constraints. Now, 
the analysis can not be made explicit in the most general situation. Answers involve certain 
summations or integrals which are well defined, but whose properties are not so transparent. 
To make explicit comparisons between properties of \l/ 7 and ^ hy , we consider two specific 
examples which have drawn a great deal of attention in connection with quantization of 
non-trivially constrained systems and 'the issue of time'. In both cases, the group averaged 
states #P hy are again semi-classical. In fact, in the more non-trivial case when the reduced 
phase space is compact, group averaging actually improves semi-classicality by reducing the 
fluctuations! These results are non-trivial in that, in the current analysis, they arise from 
delicate cancellations between terms which can not be foreseen until the very last step. It 
is these subtleties that suggest that there may well be a more general underlying structure 
yet to be unveiled. 

We will conclude this Introduction with two remarks. First, while there is considerable 
literature on the coherent state quantization [l[ 0, H 0, HI , the focus there is on constructing 
a quantum theory of non-trivial systems using coherent state techniques. Our emphasis 
is different. We are primarily interested in semi-classical issues in the resulting quantum 
theories. Secondly, our goal here is not to provide a complete and exhaustive analysis of the 
application of the group averaging technique to construct candidate semi-classical physical 
states. We only wish to present a few simple examples in the hope that the explicit and 
rather intriguing results in these cases may stimulate further research. 

II. PRELIMINARIES 

This section is divided into two parts. In the first, we fix our notation and recall certain 
facts about standard coherent states that will be used repeatedly in sections |IV| and El 
In the second, we summarize the group averaging technique which forms the basis of our 
construction. 

A. Linear phase spaces and coherent states 

Consider a system with D degrees of freedom, with linear phase space T = R 2D . It is 
convenient to use a canonical basis so that the position vector 7 of each phase space point 
can be specified by its 2D components (qi,Pi), with i — 1, 2, . . . , D. In examples, the basis 
will be adapted to the constraints at hand. The symplectic tensor Q on T is then given by 

fi(T> 7) = PiQi ~ QiPi > (!) 

where we have used the summation convention for repeated indices. (Unless stated otherwise, 
we will follow this convention throughout this paper.) 

In sections IIVI and El T will serve as the kinematical phase space. Since it has a linear 
structure, the kinematical quantization can be carried out in a standard fashion. In the 
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discussion of coherent states, it is convenient to use Fock quantization. This requires the 
introduction of a Kahler structure on T, which in turn requires new scales to enable one to 
define dimensionless holomorphic coordinates Zi. For definiteness, we will let the new scales 
£i have dimensions of length and adapt the initial choice of the canonical basis to the Kahler 
structure so that Zj is given by 

Ql +i^p- (2) 



(where there is no summation over i on the right side). The length scales £i are often 
provided by the physical parameters of the system under consideration. For a single har- 
monic oscillator, for example, £ = ^Jh/muj. The Kahler metric endows the position vectors 
7 = (qi,Pi) in T with a norm, given by 

\j\ 2 = ZiZi- (3) 

In quantum theory, Zi are promoted to annihilation operators cij with commutation rela- 
tions [aj,a]] = 5ij 1. Given a specific point a in T with coordinates (q°,p°), following (J2J) we 
will set 

«.-=^ + i^ (4) 

" V2£t V2h [) 

(Again, there is no summation over i on the right side.) Then, the canonical, normalized 
coherent state ^ a which is 'peaked' at a is given by 

= e''/ ^ \n 1 ,n 2 ,...n D ) ■ (5) 



„ /n]T • • • y/n D \ 

In the configuration representation, this state can be expressed as 

»•(«) -h e ^m^^ («> 

i=i ^ * ' 

In these states, the expectation values and uncertainties of the basic canonical variables are 
given by 

(*«| ft |*„) = q°, (*«|pi|*«> =P« (7) 

and 

(Aft) 2 = (*«| g 2 |* a ) - [(* a | ft |M/ a )] 2 = \ £\ , 

(A^) 2 ^(vI/ a |p 2 |vl/ a )-[(M/ a |p l |vI; Q )]2 = itf/^. (8) 

Thus, for each value of i, the product of uncertainties (Aft)(Apj) is minimized. In fact, 
the requirement that (jZJ) and (jHJ be satisfied by a state ^ a suffices to restrict ^f a to be a 
coherent state, given by (jnj) modulo a phase factor; this provides an alternate and more 
physical definition of a coherent state. Note also that the length scales £j we introduced 
in the specification of the Kahler structure have a direct physical interpretation in terms 
of coherent states: £i/\/2 are the uncertainties in ft in the coherent states ^f a , canonically 
defined by the Kahler structure. 



4 



Finally, for future use, we note two properties of coherent states which can be easily 
verified. First, the scalar product between two of them is given by 

(tf^a) = e-(H 2+ ^ 2 )/ 2 e? a . (9) 

Second, given any function on the phase space that can be expressed as a polynomial 
F(zi,Zj), the corresponding normal ordered operator :F(a\,aj) : has the following expecta- 
tion value: 

(^| : F(a\, a,) : |tf a ) = a,) (V^*) ■ (10) 



B. Group averaging 

For a large class of physically interesting systems, the group averaging method provides a 
technique to extract physical states starting from kinematical ones. In this sub-section we 
will present a pedagogical summary, focusing only on those aspects of the procedure that are 
central to sections ll V I and IV! In particular, we will consider a single constraint, although all 
our considerations generalize in a straightforward fashion if the set of quantum constraint 
operators is Abelian. For a discussion of more general cases and a treatment of subtleties 
and technical caveats see, e.g., Ref. 

Let us first suppose that the constraint operator C on the kinematical Hilbert space 7Ykm 
is self-adjoint, and that zero is a discrete point in its spectrum. Then, the kernel of C is a 
subspace of 7ikin- Therefore, to extract physical states, one just has to project kinematical 
states to this subspace. In the case when the 1-parameter group U (A) = e~ lAC generated 
by C on Ti^m provides a representation of U(l), the projection procedure can be explicitly 
carried out through an integration: Given any \l> G 7ikin, set 

% hy :=j£d\e- [xd y, (11) 



e~ lAC = 1. Then ^ p h y satisfies the constraint, C ■ ^ p h y = 0, 
Since these physical states belong to Tikin, the scalar product 
In terms of the general 'seed' kinematic states ^ used in ([lljh 

(^ phy |$ phy ) = (P*|P$) = = \j dA(e- iA ^|$> . (12) 

Thus we have reformulated the expressions of the projection operator and the inner product 
using averages over the group generated by the constraint. This reformulation can be readily 
carried over to the non-trivial case where the kernel of the constraint operator does not 
belong to 7ikin- 

Let us now consider the more difficult case when U(X) = e~ lXC provides a representation 
of the group K on H^ n . In this case, one starts by choosing a suitable dense sub-space S of 
^kin (specified below). Then, for each $ e 5, one defines a 'generalized bra' via 



where A is chosen such that 
and is thus a physical state, 
between them is well-defined, 
this can be expressed as 



(* phy |:=^ J d\(e- ixd V 



(13) 



5 



where the integral is now over the entire real line and K is a constant (independent of 
whose value can be chosen conveniently. While for any fixed A, {e~ lXC \?| is a well-defined bra, 
the integral fails to have a finite norm in 7Ykm- The round bracket in {^ p h y \ emphasizes the 
fact that now the physical states do not belong to H-km- Rather, they represent 'distributions' 
in the following sense: the sub-space S of Ti^n is chosen such that (^ p h y | is an element of 
the (topological) dual S* of S (with respect to a topology which is finer than that of 7ikin)- 
The meaning of (|13|) is simply that the action of this element of S* on any $ of S is given 
by 

(* phy |$> = 1 J dA(e- iA6 *|$). (14) 

Since ^ p h y no longer lies in 7ikin, one now needs a new prescription for defining the scalar 
product between two physical states. The idea is to use expression (JT2j) as a motivation and 
set the physical inner product to be 

(^phyl^phy) = (^phy|$> , (15) 

so that the norm of physical states is given by 1 

ll*p hy || 2 = ^ / dA(e-^*|*). (16) 

Thus, in the present case when the group generated by C of 7ikin is non-compact, the group 
averaging procedure provides the inner product in the physical Hilbert space only up to 
an overall multiplicative constant K. Its value is generally chosen to remove an overall 
irrelevant constant from the expression of the physical inner product. 

For a general choice of the initial subspace S, there is no guarantee that the norm (fT6*|) 
would be finite and positive. The 'art' in the group averaging procedure lies in selecting a 
dense subspace S of 7ikin such that: i) the right side of (fT%Jl is well-defined for all \I>, $ 6 S] 
i.e., typhy is a well-defined distribution over S; and, ii) the norm (jTfij) of each ^ p h y is non- 
negative, vanishing if and only if ^/ p h y vanishes. The procedure succeeds in its goal of 
constructing the physical Hilbert space only if such a S can be located. If more than one 
viable candidate exist, the resulting quantum theories could well be inequivalent. However, 
for the semi-classical issues discussed in this paper the situation is simpler, in that there is 
a natural choice which can be shown to be viable: we will let S be spanned by finite linear 
combinations of kinematical coherent states. 

For later use, we introduce the following notation: the physical expectation value of an 
observable O will be denoted by (C%h y = (^ / P h y |C > ^ / )/( v I / ph y |^')- In the explicit examples, 
we will encounter both types of constraints: in some cases U(\) = e~ lXC will provide a 
representation of U(l) and in others of R. 

Remark: Klauder []J has introduced an alternative procedure to quantize constrained 
systems. When the group generated by constraint operators on Hkm is non-compact, the 
idea is to extract physical states in two steps. In the first, one projects on the subspace of 



This procedure can be heuristically understood as follows. (|13fl extracts from VP a physical state ^ p hy S 
S* . This extractor E can be formally thought of as — 5(C)^ . Therefore, the naive definition 
(^phyl^phy) = (8{C)^\5(C)^i) of the norm that one may first think of is divergent. In the correct 
definition, one of the two delta-distributions is simply dropped, thereby removing the obvious infinity. 
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7Ykin corresponding to a portion [—5, 5} of the spectrum of the constraint operator. In the 
second, one takes the limit as 5 goes to zero. A case-by-case study is needed to check if 
the second step can be carried out successfully. This is analogous to the fact, in the group 
averaging procedure, that there is no a priori guarantee that the candidate physical inner 
product would be positive definite. However, to our knowledge, the precise relation between 
the two procedures has not been studied. 



III. SEMI-CLASSICALITY 

Even for simple unconstrained systems, the textbook treatments of semi-classicality gener- 
ally overlook two points which are important to our analysis. We will first discuss these, 
and then consider issues relevant to constrained systems. 

Fix a point a in a linear phase space F with coordinates (q°,p°)- Our task is to spell out 
what we mean by a semi-classical quantum state which is 'peaked at this classical state'. 
The intended meaning is intuitively clear and, although it is not always stated explicitly, 
one generally has the following idea in mind: a semi-classical state $ a should be such that, 
for all well-behaved functions F(qi,pi) on phase space, the expectation values (& a \F\Q a ) 
are close to F(q°,p1) and the fluctuations small. However, such semi-classical states simply 
don't exist unless the class of observables is greatly restricted. 

Let us discuss this issue in some detail. Take a harmonic oscillator. Then, one generally 
takes coherent states ^ a of (0) as representing semi-classical states peaked at the point a of 
the phase space T. Indeed, if one restricts oneself to the set, say, of three key observables, 
q, p, H of the system, where H is the Hamiltonian, ^ a would satisfy the above criteria (if 
the words 'close to' and 'small' are interpreted appropriately; see below). However, if the 
set also includes the observable e H ^ e with e < ftw, coherent states would strongly violate the 
criteria. Of course, for the harmonic oscillator the new observable is rather strange and it 
is difficult to justify its inclusion in the list on physical grounds. 

However, there is in fact an interesting situation in which the analogous observable is 
of direct physical interest. This is provided by the quantum theory of the Einstein- Rosen 
waves in 4-dimensional general relativity 0] . This system is equivalent to the axi-symmetric 
sector of 2+1 dimensional general relativity coupled to a scalar field (l0|. In the Hamilto- 
nian framework, one can arrange matters so that the true degrees of freedom are coded in 
the scalar field <f> propagating on a fiducial 3-dimensional Minkowski space. The physical 
scalar field can be identified with <p and the physical metric coefficients are also completely 
determined by it. Outside the support of 0, the only non-trivial metric component g is given 
by g = e GH , where H is the Hamiltonian of the scalar field in Minkowski space and G is 
Newton's constant in 2+1 dimensions. Since the true degree of freedom is in 0, one can carry 
out a Fock quantization and represent all interesting observables, including g, as operators 
on this Fock space. Let us now consider the problem of defining semi-classical states for 
this system. If the list of observables of interest includes only the (smeared) scalar field 
and its (free-field) Hamiltonian, the standard coherent states are good semi-classical states. 
However, if the list includes also the metric component g, they are not Conversely, one 
may construct states which are semi-classical for g but these then fail to be semi-classical 
for the (smeared) field Thus, there is no canonical notion of semi-classicality for the 



system, independent of one's choice of observables. 

The first lesson then is that to ask for semi-classical states, one must first specify a class 
of observables for which the states are to be semi-classical. A state may be semi-classical for 
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one choice but not for another. 2 

The second subtlety has to do with the notion of fluctuations. The requirement that the 
fluctuations of an observable F in a state \l/ be small is generally formulated as 

(^F)i _ (y\F^}-mP\n 2 5 2 (17) 

K^|F|^)| 2 \(^\P\^)\ 2 

where 5 denotes the 'tolerance' one wishes to allow. There is, however, the following problem 
with this proposal: if the expectation value (F)$ vanishes, the requirement can never be 
met. 3 For constrained systems this issue is especially serious because we would, in particular, 
like to consider kinematical coherent states which are peaked at a point on the constraint 
surface with only a small spread and compare them with physical semi-classical states. 
Criterion (|T7|) would forbid us from considering such states. 

The following simple example shows that, for a number of semi-classical considerations, 
(fTTj) is indeed an incorrect way to encode the idea that the fluctuations of F are small. Con- 
sider a macroscopic, 1-dimensional harmonic oscillator, such as a pendulum, with position q 
as one of the observables used to test semi-classicality of states. Consider a point (q°,p°) on 
the phase space, with real and imaginary parts of a (of Eq (@J)) large compared to 1. Then 
the coherent state ^ a of (jSJ) is a good semi-classical state by any reasonable criterion (and 
also satisfies (|17j) . Now, as is well known, under time evolution the coherent state remains 
coherent and its peak simply follows the classical trajectory of the oscillator. Furthermore, 
the uncertainty Aq is time independent. Yet, since the classical trajectory passes through 
phase space points at which q vanishes periodically, criterion ()17|) would have us say that 
the state violates semi-classicality at those times. Clearly, this is just wrong! Physically, 
since the state is semi-classical initially and does not spread, it is semi-classical at all times. 
Thus, we need to modify the criterion (|T7j) . 

Let us use direct physical considerations to develop the appropriate replacement. Given 
any state, if one is allowed to make arbitrarily accurate measurements of any one observable, 
one would invariably find deviations from the classical behavior. Thus, a quantum state can 
be well approximated by a classical one only if the experimental accuracy is limited. To 
test semi-classicality, we must supply information about these experimental limitations, i.e., 
tolerances which are fixed at the outset. We will need two sets of numbers, one specifying 
the tolerance in the accuracy of the expectation value, and the other one that in fluctuations. 

These considerations lead us to a specific notion of semi-classicality that will be used in 
sections HVI and IV! A state \l/ a will be said to be peaked at the point a G T and semi-classical 
with respect to a given set of observables Fi if 

|<tf a |i?i|tf a ) - Ft(a)\ < et and (AFj 9a < $ , (18) 



2 Heuristically, this can be understood in the following terms. The phase space structure provides a natural 
symplectic measure and, in the Bargmann-type representation, semi-classical states peaked at a point 7 
are concentrated in a neighborhood of size hP . However, the 'shape' of this neighborhood can vary. For 
one shape, they would be good semi-classical states for one set of observables while for another shape, 
they would be semi-classical for another set of observables. We thank Carlo Rovelli for this remark. 

3 One might think that this problem can be trivially overcome simply by replacing F with F + cl for a 
suitable constant c (of appropriate dimension). This strategy has two problems. First, by choosing the 
constant to be sufficiently large, any state can be made to satisfy l|17|) . More importantly, the example 
below shows that one continues to run into a problem no matter what constant one adds. 
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where and Si are pre-specified tolerances determined by the desired experimental accuracy. 

Finally, let us consider constraints. As explained in section HI the idea is to use the group 
averaging technique to extract physical semi-classical states \I/P hy — i.e., semi-classical states 
which are annihilated by the constraint operator — starting from standard coherent states 
^ a in 7ikm- Since the notion of semi-classicality is relative to a set of observables, we will 
begin by fixing the set of Dirac observables Oj of interest, together with tolerances e; and 
Si- By definition, the physical states \1/P hy will be semi-classical if they satisfy (|18|) . The 
issue then is: Can we make a suitable choice of that will guarantee that the \I/P hy are 
semi-classical? An example of a sufficient condition for the answer to be affirmative is 

1(0^- (^10,1^)1 < | ei and KAd^y-CAd^KH- (19) 

For, if this were the case, we would just need to select the kinematic coherent states \I/ a 
to satisfy our criterion (fTHj) for Oi, with tolerances e«/2 and Si/ 2 (assuming the Si/ 2 are 
compatible with the uncertainty relations). In the next two sections we will see that these 
conditions are met, even in cases where there is no obvious a priori reason for this to happen. 

IV. LINEAR CONSTRAINTS 

Let us now consider constraints of the type 

C := KiQi + K iPi - A = , (20) 

where Ki, Ki and A are any real constants. For simplicity but without loss of generality we 
will assume that C is dimensionless. We will first consider a single constraint (J2(J|) and show 
that group averaging of the kinematical coherent states ty a provide physical semi-classical 
states. This discussion can be easily extended to incorporate a set of commuting linear 
constraints. In the second part of this section, we will illustrate these constructions by 
applying them to the Gauss constraint of the quantum Maxwell field in Minkowski space. 4 
Constraints of linearized gravity can be treated in a completely analogous fashion. 

A. General structure 

Using the definition (J2J) of the complex coordinate z iy the constraint ()20|) can be written 
as C := RiZi + KiZi — A = 0, where the complex numbers ^ are related to Ki, Ki in the 
obvious manner. Given any coherent state ty a in the kinematical Hilbert space, where a is 
not necessarily on the constraint surface, it is easy to verify that 

U(\)* a := e~ ixd V a = e" iAC ^ * a(A) , (21) 

where C(a) is the value of the classical constraint C at a and <x,(A) = a,j — iXnj with <x,- 
labelling the initial phase space point a as in (jl}. 



4 However, since our primary interest lies in the relation between kinematical and dynamical semi-classical 
states, we will refrain from digressing into infinite dimensional subtleties. Readers who are familiar with 
Gaussian measures on infinite dimensional spaces should be able to fill in the details easily. Those who 
are not so familiar, can follow the reasoning using the close analogy with the finite dimensional example 
discussed in section ITY Al 
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Thus, apart from a phase factor, the image of the coherent state \l/ a is another coherent 
state ^a(A) whose peak is displaced by — i\Kj. The physical state \I/P hy is just a (continuous) 
superposition of these displaced coherent states. 

To make the properties of \I/P hy and the expressions of Dirac observables transparent, it is 
convenient to tailor the initial choice of canonical coordinates (qi,Pi) to the given constraint. 
It is obvious that we can always orient our basis so that we have Kiqi + KiPi = qi, and the 
constraint reduces then to the simple form 

C = q x - A = . (22) 

For simplicity, let us choose all ^ to be equal. Then in the g-representation, the action of 
U (A) on ty a further simplifies 

U(\)* a (q) = M a e iAA j^+P°-i/K) e-W-WM* , (23) 

where J\f a — e~ ip °' q °/ h /{ixC 2 ) / 4 is the normalization constant of the initial coherent state, 
and \q — q°\ 2 = (qi — — qf)- Thus, now the only shift in the peak of the coherent state 
is in the first component of the momentum: p\ —>■ p\ — A. Consequently, (apart from an 
overall constant phase factor) the operator U(X) simply moves the peak of the coherent state 
along the gauge orbit generated by the classical constraint in the phase space T. This direct 
interplay between classical and quantum theories is tied to the fact that the constraint is 
so simple. The physical state can also be readily calculated. We fix the overall rescaling 
freedom in ()13|) by setting K = 2ir / (ni 2 ) 1 ^ 2 for later convenience, and we obtain 

(ni 2 ) 1 / 2 r 

^(q) ■= J d\U(\)* a (q) 

= (7r£ 2 ) 1 / 2 ^ a 5(q 1 -A)e^°-^ n e^-^ 2i2 

= (n£ 2 )^ 2 5( qi -A)^ a (q). (24) 

Thus, just as one would have expected from section III Bl \I/P hy (g) is a genuine distribution; 
it fails to belong to the kinematical Hilbert space 7ikm- Nonetheless, because the constraint 
C is diagonal in the g-representation, the group averaging procedure has a simple interpre- 
tation: modulo a state- independent multiplicative constant, the physical state is just the 
restriction of the initial coherent state to the constraint surface q± — A. 

To treat more complicated situations such as the Maxwell theory discussed in section 
II V Bl it is useful to spell out the interplay between the reduced phase space quantization 
and the group averaging method used here. Since the constraint q\ — A = is so simple, there 
is a natural projection from the phase space T and the the 2(D — 1) dimensional reduced 
phase space V, spanned by q T ,Pp where / = 2,...D. Therefore, with every kinematical 
coherent state one can naturally associate a coherent state ^ in the reduced phase space 
quantization. Finite linear combinations of ^^ y , where a lies on the constraint surface, span 
a dense sub-space of 7Y p h y - For these physical states, it is straightforward to check that the 
inner products (^P hy |^P hy ) in the group averaging method are equal to the inner products 
Oi^lNj/g) in the reduced phase space quantization: 

(^ hy |< y > := = • (25) 

Furthermore, in this example, the physical inner product is also simply related to the 
kinematical inner product in the following sense. Suppose as before that a, (3 lie on the 
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constraint surface. Let a, (3 be points on a 'gauge fixed surface' (so that the corresponding 
configuration and momentum coordinates satisfy g° = A and p\ = A' for some fixed A'), 
which are related to a, (3 via olj = &/, fa = fa for I = 2, . . . D. Then, we also have: 

<^ hy K y > = (*a|^> , (26) 

where the inner product on the right side is that in the kinematical Hilbert space. Note, 
however, that this equality of inner products does not imply that are physical states; 
they do not satisfy the constraints. As emphasized above, because of their distributional 
character (see (j22)), none of the physical states belong to the kinematical Hilbert space 
while clearly do. Nonetheless, ()26|) provides a convenient tool to evaluate scalar products 
between the physical semi-classical states. In particular, it implies that the physical states 
under consideration are normalized: 

||*WH 2 = / d D q^(q)y a (q) = l. (27) 

This is a consequence of our choice K = 2n/(n£ 2 )^ of the overall constant K in f)13|) . 

Let us now turn to observables. The form (|22|) of the constraint immediately provides 
us with a complete set of Dirac observables which strongly commute with C: qj,pi, where 
/ = 2, 3, D. We will work with a larger set, consisting of general polynomials F(zj, zj) 
and their normal ordered quantum versions F = : F(a\,aj) :. As noted in section ITT) the 
expectation values of these operators in the kinematic coherent states ty a are just the values 
F(ai, otj) of the classical functions F, evaluated at the points a of the phase space: 

(* a \F\* a ) = F(a I ,a J ) . (28) 

To calculate fluctuations, by moving all annihilation operators in the expression of F to 
the right of all creation operators and keeping track of the commutator terms that result, 
one can express F as a linear combination G of products of normal ordered creation and 
annihilation operators. Hence, the fluctuations are given by 

(AF) 2 a = (F 2 ) a - ((F) a ) 2 = G(aj, aj) - (F(aj, aj)) 2 . (29) 

Since F and G do not involve a± and a\, it is easy to calculate the expectation values and 
fluctuations also in the (normalized) physical states \1/P hy (x). One obtains 

(F)^ = F(aj, on), and (AF^) 2 = G(a h aj) - (F(a 7 , aj)) 2 . (30) 

Thus, the expectation values and fluctuations of our large class of Dirac observables are 
identical in the kinematical coherent states and in the physical states obtained from them 
by group averaging. Therefore, if we ensure that the are semi-classical for the given set 
of Dirac observables, the \1/P hy will also satisfy our semi-classicality criteria. 

Given the form of \1/P hy and the form of Dirac observables, this overall result could have 
been anticipated. Nonetheless, there are two aspects which are rather surprising, at least at 
the outset. The results hold even when i) the initial choice of ty a is such that the fluctuation 
in the constraint is arbitrarily large; and, ii) the point a does not even lie on the constraint 
surface. 
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B. Example: Gauss constraint for the Maxwell field 



To make the general construction of section IIV Al concrete we will briefly discuss a field 
theory example in which there is an infinite number of commuting, linear constraints. Also, 
to provide a complementary perspective, we will use the Fock rather than the configuration 
representation. 

Let M denote a t — constant plane in Minkowski space-time. The kinematical phase space 
T of the Maxwell theory consists of pairs (A a (x), E a (x)) of suitably regular fields on M. The 
constraint surface is defined by D a E a = 0. In the Fock quantization, it is more convenient 
to work with the Fourier components of these fields. Let us introduce an orthonormal basis 
(k a , m a , rh a ), where k is the unit radial vector and m a , rh a provide a normalized 'spin-dyad' 
on each 2-sphere to which k is normal. Then, the phase space can be coordinatized by pairs 
(qi(k), Pi (k)) of fields in momentum space, given by 

A a {x) = * /2 J d 3 ke ik x (qi{k)k a + q 2 {k)m a + q z {k)m^j 

E a (x) = -j^y^ j d 3 ke ik - x ( Pl (k)k a +p 2 (k)m a + P3 (k)m a y (31) 

The Poisson brackets {E a (x), At>(y)} = 5%5 3 (x,y) then imply that the fields qi{k), Pj{k) are 
canonically conjugate in the sense that 

{q i (-k),p j (k')} = 8 ij S 3 {k,k') . (32) 

The standard Kahler structure is given by the positive and negative frequency decompo- 
sition. The holomorphic coordinates Zi(k) are now given by 

Zj (k) = {y\k\q 3 {k) - -2— Pj {k)^ . (33) 

The Gauss law D a E a = is equivalent to pi(k) = which, in turn, can be recast as an 
infinite set of commuting constraints, 



Cf(k) := J d 6 kf(k)(z 1 (k)-z 1 (k)) = Q, (34) 

one for each regular function f{k) in the momentum space (e.g., an element of the Schwartz 
space in R 3 ). Together, these constraints are equivalent to the requirement that E a have no 
longitudinal modes P\{k). A complete set of Dirac observables is therefore given by arbitrary 
real-valued functions of the transverse modes zi(k), zi(k), where I — 2, 3. 

The creation and annihilation operators correspond to the classical functions Zi(k) and 
Zj(k) respectively. They satisfy the usual commutation relations: 

[d t (k),a](k')] = 6 l:i 6 3 (k,k')i. (35) 

The kinematic Hilbert space H^n is the Fock space obtained by operating repeatedly with 
the creation operators on the vacuum state |0). For each point a in the phase space, a 
coherent state ^/ a , peaked at Zi{k) = a^k), can now be constructed in Hkm following the 
procedure outlined in section ILTAl 

1^ _ e /(d s fc/|fe|)(a(fe)-at(fe)-a(fc).a(fe)) ^ _ (3^ 
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For each mode k, we now have a linear constraint pi(k) = 0. Hence, the passage to 
quantum theory is structurally similar to that in section II V Al Let us therefore begin 
with a kinematic coherent state ty Q and apply the group averaging procedure to it. The 
physical state (\I/g hy | is a well-defined distribution over the sub-space S of the Fock space 
7Ykin spanned by finite linear combinations of coherent states. It does not belong to 7Ykm 
because its kinematical norm diverges. However it is straightforward to calculate the action 
of the distribution (^Q hy | on elements |\Jjg). In the interesting case when a and j3 lie on the 
constraint surface, the action is simple to write down (see (j2SI)) : 

(tffW = (*ai^> , (37) 

where ai(k) = 0, 6ci(k) = ai(k) and $i(k) = 0, $i(k) = (3i(k). This equality is useful in 
computations: ^^ y , with a lying on f span a dense subspace of Tt p h y and the physical 
inner product between these states is given by: 

(* p « y K y ) = (n\*0) ■ (38) 

However, as emphasized in section II V Bl after (|26[). it does not imply that physical states 
form a sub-space of 7ikin- While states in Hkin are peaked at points on the constraint 
surface of the phase space, they have fluctuations also away from the surface. They do not 
satisfy the quantum constraint. 

Since Dirac observables O in our complete set depend only on zj, it follows from the 
arguments given in section II V Al that the expectation values and fluctuations of O in 
and \I/P hy are identical. 



V. QUADRATIC CONSTRAINTS 

We will now consider constraints for which f is the level surface of a quadratic function on 
T. As in section HTH without loss of generality, we will suppose that C is dimensionless. Fur- 
thermore, for technical simplicity we will make a further restriction: the constraint function 
will be assumed to be of the form: 

C{qu Pi) ■= Sij qiqj + A Sy + A i4 q^j - A = , (39) 

where is a symmetric matrix, Aij an anti-symmetric matrix, A a constant with dimensions 
[L 2 / (Action)] 2 , and A a real constant. This class includes a number of interesting cases. For 
example, in geometrodynamics, the function obtained by smearing the diffeomorphism 
constraint with any vector field on the 'spatial' manifold M, 

Cjf(q,p)= [ P ab (x)£^q ab (x)d 3 x , (40) 

J M 

is of this type. (In the 'super-index' notation introduced by Bryce DeWitt, i = (x;a,b), 
this function is of the type = qiNijPj where Nij is an anti-symmetric matrix. Thus, in 
this example, = and = Nij.) As we will see in sections I V CI and IV Dl some of the 
finite-dimensional examples that have been studied extensively in the literature are also of 
this type. For simplicity, in the detailed analysis we will consider a single constraint, but 
it is rather straightforward to extend it to allow for a set of commuting constraints each of 
which is of type (JBTJj) . 
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A. Setup and physical states 



), the constraints we consider can be 



In terms of the holomorphic coordinates 2, of Eq 
written as 

C(qi,pi) := KijZiZj - A = , (41) 

where Hermitian matrix (and we have chosen all li equal to one another in the 

definition |21 of holomorphic functions Zj). As in section II V A[ the analysis becomes more 
transparent if the initial choice of canonical coordinates is adapted to the constraint at 
hand. Let us therefore choose (qi,Pi) such that the Hermitian matrix Kij is diagonal, with 
eigenvalues Ki (which are of course real). Then, using as before normal ordering, the quantum 
constraint operator becomes 

C = Kj Nj - A i , (42) 

with Nj the jth number operator, Nj = fit a, (where there is no summation over j). 

The action of U(X) '■= e~ lXC on the kinematical coherent states ty a can be calculated in 
a straightforward manner: 



-i\C 



D 



iAA —WYijKja-aj 

e e 



iAA 



D 



i=l 



i=l 

oo 



"l«i| 2 /2 



E 

rii=0 



l°i| 2 /2 



E 

rii=0 



"ru 



[OLi 



"nv. 



\ni) 



e UA |*a(A)), 



(43) 



with atj(X) = e 



-iAft. 



a, 



(where there is no summation over j). Thus, apart from a constant 



phase factor, the image of $f a under U(X) is again a coherent state, the peak being shifted 
from a to a(X). It is easy to verify that, on the classical phase space T, a(X) is precisely the 
gauge orbit passing through a, generated by the constraint function C. Consequently, as 
in section IIV[ there is a close interplay between classical and quantum theories also for the 
quadratic constraints now under consideration: the action of U(X) simply moves the peak 
of the kinematical coherent state along the gauge orbit generated by C on F. 

Since we wish to use the group averaging technique, we need to restrict ourselves to the 
case in which U(X) provides a representation of a group on Hkin- The expression (f4l?J) for the 
action of U(X) shows that it provides a representation of U(l) if and only if there is a real 



number A such that e 



iAA 



-iA« 



1 for all i. We will consider constraints of this type. 



This in turn guarantees that the kernel of C, the physical Hilbert space 7i p h y , is a subspace 
of T^kim and group averaging now falls in the simpler of the two cases considered in section 
111 Bl defining a projection operator P from 7Ykin to H p h y . The condition also means that 
all ratios k«/A and have rational values, and since we can multiply C by a constant, 

we may take all Ki as well as A to be integers, and A = 2n (choosing the integers with the 
smallest absolute values fixes them uniquely up to an overall sign). In addition, it can be 
seen from the form (}4*2*|) of C that, in order for C to have a non-trivial kernel, Ki and A 
must be such that Kifii — A = for some choice of integers m, . . .no- To ensure that the 
example is interesting, we will assume this to be the case. 
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The physical state corresponding to ^ a is then given by Eq (fTTj) . or 



^ = i- [ d\U(\)* a , (44) 
Z7r Jo 

and using (|T2*|) and we find that its norm is 

' dAe -iAA e E,|«,|V A «i i (45) 



2tt 



o 



Since the integrand is a smooth function of A and the integral is over a closed interval, it 
is clearly well-defined. By expanding the coherent states in the first equation of (|4~5j) in the 
occupation number basis, the norm can also be expressed as the sum 

\2m |_, \2n D 

||*?T = e-W a £ ' ' n ,;;i <W, (46) 

ni, ...,n D 

where the Kronecker delta limits the contribution to a finite set of rii, those satisfying 
KiUi = A. In this series form, positivity of the norm is manifest, whence the result of group 



averaging is indeed a physical state. More generally, as stated in section Hi B| one can verify 
that the space S spanned by finite linear combinations of coherent states is an admissible 
dense sub-space of 7Ykin to serve as the 'seed' in the group averaging procedure. 

Notice that, if the Ki are all non- negative, one can also obtain a closed form for the norm 
directly from the last integral in ([45)1 . Setting ( = e lA , the integral reduces to one over the 
unit circle in the complex ( plane, which can be evaluated using the method of residues: 



\\K hy \\ 2 = ^r- <fi dC 



«l 2 f e £>;l 2 0- e -M 2 d A 



2vri J K]=1 * C A+1 A! dC z 



e 



(47) 

C=o 



from which the result (J46|) can be recovered by evaluating the derivative. If all the Ki are 
non-positive, one can use the same procedure by initially setting £ = e~ lA . In the general 
case, however, this method gives a result expressed only as an infinite series. 



B. Physical observables 

The calculations of expectation values and fluctuations of Dirac observables can also be 
carried out in a rather straightforward manner, because the expressions involve integrals of 
the type J dA a ^x)\0\^ a ) in which both the bra and the ket are coherent states. However, 
in the general case the result can again only be expressed as a well-defined integral or a 
convergent infinite series. Properties of expectation values and fluctuations, therefore, are 
not always transparent. But in the case of linear or quadratic observables, some interesting 
results can be readily obtained. 

A linear observable is one of the form O = FiZi + FiZi, and calculations of expectation 
values and fluctuations with it are greatly simplified by the requirement that it be a strong 
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Dirac observable. From the Poisson bracket between O and the constraint (|41|) we see that 
in this case Fi — for all i such that ^ 0. In a group averaged coherent state, 



^ JO 

= J_ / dAe- iAA [F J a J (A) + F,a,](* Q(A) |^). (48) 
^ Jo 

But ^ Oj(A) = e 1KiA Oj = Oj, since for each i either F, or vanish (or both). Therefore, 

~ , /vi>p h y| |vi/P h y\ 

y = { Q ^J ||2 Q ; = 0(a) , (49) 

the classical value. Similarly, for the fluctuation we calculate 

(^y\d 2 \^) = -L /^dAe- iAA (vI/ Q(A) |(F4 + FA) 2 |^) 

JO 

= _L f d\ e" iAA [FjFj ai(A)a,(A) + 2 i^i? 04(A) ^ + 
^ Jo 

+ FiFj a iaj + FiFi] (* a(A )|*a> 
= [(Fa, + Fa,) 2 + FF] ll^f , (50) 

where the last term in the second expression arises from the commutator [a*, aj], and 
Fidi(X) = Fidti for each i, for the same reason as in the expectation value. From this, 
it follows immediately that 

(AO)l hy = (AO)L = F^ . (51) 

In the case of a quadratic observable of the general form O = ZiZi, one can use 
similar reasoning to arrive at a more limited statement. Since the Poisson bracket of O 
with the constraint (|41j) vanishes, the matrices «y- and 0^- commute. Therefore, they can be 
simultaneously diagonalized and we can write the normal ordered operator for the quadratic 
observable as O = faNi. Then 

<*pfv|0|*P*) = -]T / dAe- iAA (vI/ Q(A) |0 l aia,|^) 

= ^ V. / "dAe- iAA 0ie iAK 'a 4 a 4 (^ (A) |^ Q ) . (52) 

Z7T » Jo 

In this expression, the A dependence of the integrand cannot be simplified as in the cor- 
responding one for a linear observable, (|48|). and the expectation value (0)^ y is typically 
different from 0(a). We will analyze the difference in detail in two examples. Here, we will 
only present a general argument to show that, when the phase space point a at which \I/ a 
is peaked is far from the origin, the expectation values 

- , /vi>P h y| |vi/P h y\ 
(0)t* = [ \ ' a ; (53) 
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are necessarily close their classical values. For simplicity, let us explicitly consider the case 
in which one of the \®.i\ 2 becomes very large; without loss of generality, we may assume it is 



|a D | 2 . 



Recalling that (^f a (x)\^a) — e ' a ' 2 e^^ 2 exp ( lXK: >\ we can write the numerator of (|53j) as 



p — |a|- /"27r 

^Phyj ^ |^ph y) = e y ^ , , 2 / dA e _iA(A- Ki ) ^KIV (54) 

2tt '« y 

Each integral in this last expression is of the form (|45|) with A replaced by A — K i} and can 
therefore be written immediately in the form of the result (J46j) ; thus, 

{^\d\^)=e-^j^L\a^ £ |ai|2ra V- > f ng ^,A- Kl V (55) 

^ \ ^ yii!. ..Tin'. I 

i=l \ ni,...,no ' / 

Both the numerator (jo3j) and the denominator (J4T?)) of the expectation value are finite poly- 
nomials in the |aj| 2 , which implies that when \old\ 2 — * oo the dominant terms in the sum- 
mations over m, ■■■,ri£) are those with the highest powers of |a_o| 2 - Thus, to determine the 
asymptotic value of (O)^ 7 we need to find the highest values of no contributing to the sums 
in (gUl and (gB| . 

Call Nn the highest value of rip contributing to the norm ()46|) . The classical constraint 
implies that, as \ac\ 2 — > oo, A ~ K D \an\ 2 ] together with the restriction on the n i; this then 
implies that Nn ~ |«d| 2 - Now turn to (|55|) and consider first a term with i ^ D. Depending 
on the values of all the Kj, the highest value of rin in it may be either or a smaller number; 
in the former case, the whole term is of order 1 in |«_d| 2 , in the latter it has a decreasing 
behavior. For the i = D term in ()55|) the situation is simpler, since the highest value of no 
in it is always Nd — 1; in fact, the leading term in the sum is |a£>| 2< - Ar - D ~ 1 ' ) /{Nd — 1)' times 
the same combination of the remaining atj and rij that multiplies |a£>| 2ArD /iV£>! in the norm 
(US)). Thus, as |«d| 2 — > oo, 

,^ phy (^\d\n hy ) , , |2 \a D \^ ND - l) l{N D - 1)! 2 

Two examples discussed in detail below will illustrate this general feature of our physical 
states. For some quadratic observables ((9)£ hy will actually coincide with the classical value 
0(a). These examples will also show that the fluctuations of quadratic observables in the 
physical states differ from the corresponding fluctuations in the kinematical states even for 
the simplest observables, but do so in a controlled fashion. Interestingly, in some cases they 
are smaller than the latter. 



C. Example 1: Constrained total energy 

In our first example, we will consider two coupled harmonic oscillators of the same mass and 
spring constant, subject to the constraint that the sum of their energies have a fixed value: 

C-E.jl + ^-^O. (57) 
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This example has drawn attention in the literature (see, e.g., Refs (HI 0, H|) because a 
constraint that fixes the total energy provides a toy model for the Hamiltonian constraint of 
general relativity in the spatially compact case (where the energy is zero), and can be used 
to study the associated 'problem of time'. 

To conform to our general framework, we need to divide C by a constant with dimensions 
of energy to obtain a dimensionless constraint function. While any constant will do, we will 
use the natural choice hu. Then, using for I in the expression (J2J) of z\ the natural length 
scale t = ^yh/muj, the constraint can be re-expressed as 

C := — C = Zl zi + z 2 z 2 - A = , (58) 

where A = A/hw. The kinematic phase space T is IR 4 ; the constraint surface f is a 3-sphere; 
and the gauge orbits generated by the constraint function C provide a Hopf fibration of F. 
Thus, the reduced phase space T is a 2-sphere. Because of the topological non-triviality, 
although T is 2-dimensional, we need a set of at least three Dirac observables to separate 
points of T. A convenient choice is 

Li = Re ziijj, L 2 = Im z\Z%, — \ {z\Z\ - z 2 z 2 ) ■ (59) 

(As the notation suggests, they are also closed under Poisson brackets, providing a repre- 
sentation of the standard angular momentum algebra. However, this property will not play 
a role in our semi-classical considerations.) 

In contrast to section HVl now the structure of Dirac observables is quite complicated and, 
as seen in the general discussion, there is no a priori reason to suppose that expectation values 
and fluctuations of these observables in the kinematic and physical quantum states will be 
simply related. To explore the relation, let us begin by constructing physical states. Our 
constraint operator is 

C = Oitt} - A = V" Ni - A . (60) 
Thus, Ki = 1 for i = 1,2. We can now write down the physical coherent states 

-(|ai| a +|a a | a )/2 rto °° n n m 

= ~ — / dA £ r^rb 6 ■ ^ 



o — ^ n yn\ym\ 

n,m=0 v 5 



Using the fact that the Fock basis \n,m) is an eigenbasis for the constraint, satisfying 
C\n, m) = n + m — A, where A = k := n + m is an integer, we get 

l*f y ) = 7^ / 2W dAe iAA |(« ie - iA ), (« 2 e" iA )) . (62) 
^ Jo 

Hence, from the general discussion of section IV Al we conclude that the norm of the 
physical states \I/P hy , obtained by group averaging the coherent states is given by 

(63) 
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If we now use the fact that the coherent state is chosen to be peaked at a point of the 
constrained surface f, then A = \a\ 2 = k, so we have ||^a hy || 2 = e~ k k h /k\ . Note that in 
this case, the physical Hilbert space is finite-dimensional (due to the compactness of the 
reduced phase space), with dim(7Y p h y ) = k + 1, and that the extractor operator E = P is in 
fact a projection operator on the kinematical Hilbert space, whose action is to project the 
kinematical coherent state to the subspace spanned by kets of the form \n,k — n) for a fixed 
value of k. 

We can now calculate the expectation values of the constraint and the observables Lj, I — 
1, 2, 3. 5 Let us begin with the term a\ai, and compute its expectation value (\& a | a{ai|^/ Q ) p hy 
In this case we have 

1 r 2n 

(^ a \a\ ai \^ a ) phy = — J dX e _iAA ((aie _iA ), (a 2 e~ iA ) | a\ a x a 2 ) 

= h£ /' 2W dAe -iA(A-l) e H a P e c-| Q | 2) (65) 

2tt Jo 



which can be rewritten as 



(^> a | a\ ai ^ a ) phy =\ ai \ 2 —(b ^ = \ ai \ 2 —— k . (66) 

2m J Kl=1 C {k - 1)1 



Thus, when we compute the expectation value we get 



(^oj_QlQj_^aW |„ 1 2 e fc fc-1 fc! fc 2 

; ; : — Oh — rr /C r fC — ttl 



(67) 



It is clear that the expectation value of the operator a\ a 2 will also equal its classical value. 
Finally, it is easy to show that the expectation value of the observables Lj in the physical 
coherent state labelled by (ai,a 2 ) is given, for all /, by 

(L/)phy - (Wi " 4l RFTRF • (68) 

However, since the point a lies on the constraint surface, we have | cti | 2 + |a 2 | 2 = k, whence 
the expectation values always coincide with the classical values. 

The next step is to look at the fluctuations of Dirac observables. We shall first compute 
them in the physical coherent states, and later compare them to the fluctuations in the 
kinematical coherent states. It is a straightforward calculation to show that, on a physical 
coherent state, 

<^w = ^L ( iXw + ?- (69) 



5 The expectation value of the constraint in the kinematical coherent state vanishes, and its fluctuation is 
given by 

(AC-)L = M 2 = |ai| 2 + |a2| 2 = A. (64) 

Since (AC*)£ in is proportional to E c i/hu), for coherent states with a large value of E c i, the fluctuation of 
the constraint will also be large. 
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Therefore, 



(AL, &J == (L ?)phy - «£, W )» = Lj| d ( (i^p^iy ) ^, + | . (70) 

On the other hand, the fluctuations in the kinematical coherent states are given by 

(Ai^)kin = | (|ai| 2 + \a 2 \ 2 ) = j . (71) 
Thus, the difference between the fluctuations is given by 

(AL/)L - (AL 7 g hy = i L?| d < (AC)L . (72) 

This implies that: i) The difference in the fluctuations is smaller than the fluctuation of 
the constraint operator on 7ikm; h) Group averaging actually reduces the dispersions. To 
summarize, then, if we begin with semi-classical kinematic states peaked at points on the 
constraint surface, physical states resulting from group averaging are guaranteed to be semi- 
classical. Furthermore, the kinematical calculation provides a good upper bound on the 
dispersion in the physical states. 
Let us conclude with two remarks. 

1. While the explicit calculations led to a desirable result, the 'mechanism' behind its success 
is rather obscure, even in retrospect. Indeed, from the expressions of the individual integrals 
or summations defining the numerators (\I/P hy | Lj | \I/P hy ) and the denominators (^P hy |\I/P hy ) 
in the expression of the expectation values, it is far from clear that the ratio would agree so 
well with the result on the kinematical Hilbert space. It is only after the two integrals are 
evaluated explicitly that the required cancellations occur. Thus, the desired result emerges 
somewhat surprisingly and only in the very final step. 

2. We could have chosen the quantum constraint C not to be normal ordered. Had we 
used the standard text-book ordering, even in the kinematical Hilbert space the expectation 
value of the constraint would not have been zero if the coherent states were chosen to be 
peaked at points a on f : Due to the zero point energy, we would then have ((7) kin = 1- 
Furthermore, the expectation values of the Dirac observables on group averaged coherent 
states would not have coincided with the classical values, but would have the form (Lj) p h y = 
L-i\ c \ (1 — fiu/E c i). Thus, only for large values of E c \ would the physical coherent states have 
approximated the classical values of the observables. These 'discrepancies' can be avoided 
by choosing kinematical coherent states which are peaked slightly away from the constraint 
surface, but that would have made the construction a bit obscure from a physical standpoint. 
Fluctuations would have slightly different values but share the same qualitative behavior as 
in the normal ordered case. 



D. Example 2: Constrained energy difference 

In our second example, we will consider two coupled harmonic oscillators of the same mass 
and spring constant, subject to the constraint that the difference of their energies have a 
fixed value: 
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This example has also drawn attention in the literature (see, e.g., Refs 0, |||) because it 
models some cosmological situations The system is conceptually simpler, even though 
it is more complicated from the technical viewpoint, since now the reduced phase space is 
non-compact and isomorphic to an open set in R 2 . We again have two harmonic oscillators 
with the same frequency, but now the constraint can be written as 

C= \ Zl \ 2 - \z 2 \ 2 - A = , (74) 

which is of the form (j41|) with = diag(l, —1), where as before A = A/hw. The quantum 
constraint operator has the form 

C = (a\di — — A . (75) 

(Note that, in this example, the ordering ambiguity discussed at the end of the last sub- 
section does not matter because the zero-point energies of the oscillators cancel each other 
out.) It is straightforward to write the physical coherent states, 

|v|)phy ) = ± e -(M 2 +M 2 )/2 / dA y "i tt2 e - iAC |n,m). (76) 
' 2tt J ^ o ^T^T 

Using the fact that the Fock basis \n,m) is an eigenbasis for the constraint operator, satis- 
fying C\n, m) = [{n — m) — A] \n, m), where A = k := n — m is an integer, we get 

\9**) = ± / dAe iAA e -( K |W)/2V-lE a ^ \ n ,m), (77) 

27r Jo ^ vn!m! 

but since the last term in the integral is precisely a coherent state, the physical coherent 
state can be written as 

\n hy ) = ^ /^ dAeiAA K«i e " iA )' • (78) 

Although the physical Hilbert space is now infinite-dimensional, since zero is a discrete 
point in the spectrum of C, the extractor E is again a projection operator; E = P on the 
kinematical Hilbert space Hkin- Its action is to restrict contribution only to kets of the form 
\m, k + m), for a fixed value of k. It is now straightforward to compute the norm of the 
physical coherent states, 

1 f 2n 

ll*? y || 2 := (*Jy|* Q ) = / dAe- iAA (* a(A) |¥ a ) 

Z7r Jo 

= — f 7r dAe -iAA e -(| ai |2 + |a 2 |2) e (e^|a 1 |2+e-^|a 2 |2) _ ^ 

To evaluate the integral, we can again define a new variable £ = e lA and express the integral 
over A as a contour integral over the unit circle. 

|^phy||2 = / ^-a^C+I^C 1 ). (80) 

27T1 /| f | =1 C 
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Since a is chosen to lie on the constraint surface, from the expression of the eigenvalues of 
C, it follows that A is again an integer; A = k. But now the function to be integrated has a 
pole of infinite order at the origin, whence we can not compute the integral as easily in the 
previous example. But we can still express the result in terms of special functions. Let us 
expand the integral as an infinite series 



|^rphy||2 _ e -(M 2 + M 2 ) 



a i 



m=0 



and use the identity 



£ 

n=Q 



(x/2) 



2n 



n\ (k + n)\ 



k\ (k + m)\_ 



-] ]/,■(• r) 

x 



\ai\ 



1 2m I 



«2 



1 2m 



where l m is a modified Bessel function. Then, the norm can be expressed as: 



I \[/P h y I 



-(l ai l 2+ l a2 l 2 )(|a 1 |/|a 2 |) fc I fe (2|a 1 ||a 2 |) 



(81) 



(82) 



(83) 



As one might expect from the last example, we again have three quadratic Dirac observ- 
ables: 



K Z±Zi + Z 2 Z 2 ) 



J+ '■— Z\Z2 
1 , 



Zl z 2 



(84) 



with their corresponding combinations, Ji = ^(J+ + J J) and J 2 = ^(J+ — J-)- These 
observables provide a realization of the sl(2,R) = su(l,l) Lie algebra that has been much 
studied. While this structure is completely analogous to that in the previous example, now 
there is a key difference. The reduced phase space T is no longer compact but topologically 
M 2 . One can verify that J% and J 2 suffice to separate its points. Therefore, we will first 
focus just on these observables and, for completeness, also discuss J3 at the end. 

It is straightforward to verify that the expectation values of J± — and therefore of J\,%— 
in the physical coherent states coincide with the classical values: 



(Jit) 



phy 



2 cl 



For fluctuations we obtain 



(AJ; 



ai cm 



l,2jphy 



I fc -i(2 
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\at 2 \) + Ife+i(2|ai||a 2 |) 


Ifc(2 


«i| 
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and 



Therefore, 



(aA, 2 )L 



|(|a 1 | 2 + |« 2 | 2 + 2). 



(85) 

(86) 
(87) 



(AJ; 



l,2jphy 



Ifc-i(2|ai| 


|a 2 |) + Ifc+i(2|ai||a 2 |) 


I fc (2 


ai| 


l^l) 



(88) 

In Fig. Q the fluctuations of Ji j2 in the kinematical and physical states are compared. In 
particular, the quotient 



(AJ : 



1, 2 ) phy 
(AJ; 
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l,2jkin 
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FIG. 1: The quotient [(AJ 1)2 )^ y - (AJi )2 )jy /(AJi )2 )£ in , plotted as a function of x = \a 2 \. The 
constant that fixes the energy difference is set to k = 10. Note that the ratio approaches 1/2, and 
is very close to 1/2 even when |a 2 | is not very large. 



is plotted as a function of \a 2 \ for k = 10; the qualitative behavior is the same for other 
values of k. From the numerical investigations it is clear that the quotient approaches 1/2 
very fast as \a 2 \ grows. This means that the fluctuations in both types of coherent states are 
of the same order. However, in this example, the fluctuations are smaller in the kinematical 
coherent states than in the physical states and the difference remains bounded away from 
zero. However, it is clear that if the initial kinematical coherent states are chosen with 
tolerances 2<5j/3, e,, we would be guaranteed that the group averaged physical states will be 
semi-classical with desired tolerances 5i and £j. 

Finally, for completeness, let us consider the quantum observable J 3 = ~ (a\ai + a\a 2 ). 
For the expectation value, we have: 



(Js) 



phy 



a 2 
a i 



Ifc-i(2|cti||a 2 |) , | ,2 

+ «2 



|ai/a 2 | fe I fc (2|ai||a 2 |) 



a 2 



I fc+ i(2|ai||a 2 |) 
|ai/a2| fe Ifc(2|ai||a 2 |) 



(89) 



Since |«i| 2 = k + |a 2 | 2 , one can evaluate the last expression as a function of k and \a 2 \ and 
compare it with the classical value. The quotient ({^3)ph y ) / (^3|ci) is plotted in Fig. EJ As 
can be seen, already for small values of k and |a 2 |, this quantity is very close to one; physical 
coherent states approximate very well the classical values of the Dirac observables we have 
chosen. 

The fluctuation of J3 can be obtained by considering the expectation value 



W3) phy 



\ol\ \\a 2 \ 



Ifc-i + I/c+i 



I \2i |2 I Ifc-2 + 2Ifc + Ife+ 2 

Qfi \a 2 \ 



(90) 



With this expression we can again analyze the behavior of (AJ3)p hy and compare it with 
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FIG. 2: The quotient ( JzjphyK. Js)ch plotted as a function of x = |ck2 | - The constant that fixes the 
energy difference is set to k = 10. Notice that the ratio is very close to one even for values of |o:2 1 
that are not that large. 

the kinematical fluctuation. This comparison is shown in Fig. El where the quotient 

(Ai 3 )L 

is plotted. Again the two fluctuations are of the same order. The main difference with 
the fluctuations of Ji, J2 is that now the fluctuations in the physical states are smaller but 
they very rapidly approach those in the kinematical coherent states. Therefore the physical 
states are guaranteed to be semi-classical with respect to J3. Thus, although the detailed 
behavior of J3 is rather different from the one of J\ and J2, the physical states are again 
semi-classical also with respect to J3. 

To summarize, by restricting the initial kinematical coherent states to have suitably small 
tolerances, the group averaged physical states can be guaranteed to be semi-classical for any 
specified choice of tolerances. 

VI. DISCUSSION AND OUTLOOK 

Let us begin with a brief summary of results. In section ILTTl we clarified two issues concerning 
the notion of semi-classical states. The clarifications in turn led us to a criterion under 
which states \I/P hy obtained by group averaging suitable kinematical coherent states can 
be regarded as semi-classical in 7i p h y . In sections [TV] and |V| we saw that the criterion is 
satisfied if the constraint sub-manifold T is the level surface of a linear function on T or a 
quadratic function satisfying certain conditions. Thus, the group averaging procedure offers 
a concrete and potentially powerful strategy to construct physical semi-classical states for a 
class of constrained systems. 

24 




1- 

0.8- 












0.6- 












0.4 - 












0.2- 















2 


4 


6 
X 


3 10 


12 



FIG. 3: The quotient [(AJ 3 )| in - (A J 3 )^ hy ]/ (A J 3 )£ in , plotted as a function of x = \a 2 \. The 
constant that fixes the energy difference is set to k = 10. Note that the physical fluctuations are 
smaller and very rapidly approach the kinematical ones. 

In the examples with quadratic constraints the result could not have been foreseen on 
general grounds. Indeed, even in retrospect we do not have a general understanding of 
'why' the strategy works. It is particularly surprising that, in certain cases, the group 
averaging procedure even reduces the fluctuations of Dirac observables. Now, all examples 
we considered have the feature that U(\)ty a is again a coherent state, peaked at a point 
a(X) of the phase space T on the orbit of the gauge transformation generated by the classical 
constraint C. This will not hold generally. Is this feature perhaps the key to the 'nicer than 
expected' behavior of the group averaged coherent states? For example, since the value of 
any classical Dirac observable O is constant along a gauge orbit, this feature implies that 
the expectation values (U(X)^ a \ O | U(X)^ a ) are all equal to the values 0(a) of O at the 
classical point a G T. Note, however, that it is not these expectation values that dictate the 
calculation of (\I/P hy | O | \I/P hy ). As we saw in section |VJ the calculation is governed, rather, 
by a delicate interplay between the cross matrix elements (^> a \ O \ U(X)^ a ) and the norm of 
the physical state \I/P hy . Neither of these by itself has any simple relation to the value 0(a) 
of the classical Dirac observable. Thus, while the specific property now under consideration 
of our class of constraints did simplify the detailed calculations, it does not seem to suffice to 
ensure semi-classicality of \I/P hy . Indeed, there are examples of quadratic constraints — such 
as those in the Bianchi I model — which do not share this property but where the group 
averaging procedure is again useful in constructing physical semi-classical states 14] . It 
would be extremely useful to understand the underlying mechanism which makes the group 
averaging strategy successful in a rather diverse class of examples and prove a general result 
which guarantees success for a wide class of constraints. 

The following heuristics suggest that this may well be possible. In the case of a generic 
constraint, one can simply choose the constraint function itself as the first canonical coordi- 
nate, say qi, on the linear phase space T, and then supplement it with other suitably chosen 
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functions (tailored to one's choice of Dirac observables) to form a canonical coordinate sys- 
tem (qi,Pi). Then, as in section HVl physical states would have a simple distributional 
form in the g-representation: \l/ phy (g) = 5(q , i)/((?2, <?£>)• However, in the general case, 
this canonical chart would not be adapted to the linearity of the phase space, whence it 
would be difficult to identify coherent states in this representation. Nonetheless, it should 
be possible to introduce some notion of kinematical semi-classical states (?). Then the 
group averaging procedure would lead to states ^^(q) = <5(gi)\l/^ m (g2, These are 

natural candidates for physical semi-classical states. For, Dirac observables would act only 
on the variables q2, ■ ■ ■ ,qD, whence there could be a simple relation between their expec- 
tation values in ^^ n {q) and ^^ y {q) needed to establish their semi-classicality. A number 
of non-trivial technical problems need to be overcome to determine whether these heuristic 
ideas can be made precise. Nonetheless, they indicate that there may well be a much more 
general underlying structure responsible for the success of the group averaging method in 
the few examples discussed in this paper. 
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